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NOTATION 


Fnglish Letters 


A Angstrom unit 
a equilibrium speed of sound 
B. partition function of argon atom 
B n partition function of argon ion 
A 
ce speed of light 
D diameter of cylindrical leading edge of model (aye inch) 
e, specifie lonization energy 
e, internal energy per unit mass due to translation 
h specific enthalpy 
hp Planck's constant 
K Gladstone-Dale constant 
k Boltzmann's constant 
L test section width 
M Mach number 
M. incident shock Mach number 
mM, mass of argon atom 
m mass of argon ion 
+ 
A 
mM mass of electron 
Ny number density of atoms 
N n number density of ions 
A 
N, number density of electrons 
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NOTATION (Continued ) 


index of refraction 

retracuivity 

pressure 

electronic charge 

gas constant per unit mass 

average R for the mixture 

gas constant per mole 

fringe shift; also fringe shift at A = 4500A 
fringe shift at A' = 4030A 

temperature 

fluid velocity relative to the standing shock (lab) 
atomic weight of argon atom 

average atomic weight of the mixture 


degree of ionization 


Greek Letters 


ea 


polarizability 

ratio of specific heats 

change in ( ) 

shock stand-off distance 

argon ionization temperature 

observing wave length; also A = 4500A 
4030A 


reduced mass of ion-electron pair 
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NOTATION (Continued) 


fe) mass density 

io) mean density from region 3 to stagnation line 

> argon ionization potential 

W frequency of observing light, corresponding to xX 
a plasma frequency 

Subscripts 

ak region 1, initial downstream conditions 

2 region 2, conditions entering the standing shock 
5 region 3, conditions downstream of the standing shock 
S stagnation line conditions 

A atoms 

AY ions 

e electrons 
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I. INTRODUCTION 


Treatment of high temperature shock layer flows has been of interest 
from the point of view of hypersonic and re-entry heat transfer. Fay 
and Riddell bole? developed a correlation formula for convective stag- 
nation point heat transfer rate in a dissociating diatomic gas, in 
conjunction with an experimental study by Rose and Stark [2]. Several 
investigators have applied a modification of this formula to stagnation 
point heat transfer in partially ionized argon. Reilly [3], operating 
in a regime where radiative heat transfer was considered negligible, 
measured convective heat transfer rates slightly higher than predicted 
by the Fay-Riddell equation applied to an ideal gas. Rutowski [4] also 
determined that convective heat transfer was only moderately increased 
by real gas effects in the partially lonized argon. Rutowski and Ber- 
shader [5] found that, for a shock Mach number of twelve, radiative and 
convective heat transfer rates were about equal, while the convective 
rate alone agreed with the Fay-Riddell ideal gas calculation. They were 
also able to correlate the measured radiative heat transfer to the cal- 
culated electron density, using the Kramers classical theory to relate 


radiation to the square of electron density. 


The importance of electrons in both radiative and convective ther- 
mal transport leads to interest in electron density measurement. The 
present study deals with a measurement of electron density near the 
stagnation line of a two-dimensional shock-heated argon flow around a 
cylinder-plate (Fig. 1). The purpose of the work was to make a quanti- 
tative comparison between measured results obtained by optical inter- 
ferometry and the predictions of equilibrium flow theory, using the 
Rankine-Hugoniot and Saha equations. A noble gas is desirable as the 
working fluid, since ionization can be studied without the complications 
of dissociation. In particular, argon was chosen because of its moderate 
ionization potential and relatively high molecular weight, which enabled 
significant ionization to be obtained behind the standing shock. 


*) 
Numbers in brackets subsequently designate References collected 
at the end of this paper. 





The technique of optical interferometry depends upon changes in 
index of refraction of the test gas, and has long been used to map 
densities in non-ionized flows. But Alpher and White [6] demonstrated 
that interferometry at optical wave lengths could be used to measure 
electron densities in a plasma. Since the refractivity of atoms is 
relatively insensitive to wave length, while the refractivity of elec- 
trons is highly dispersive, measurement of fringe shifts at two wave 
lengths enables the electron contribution to be separated. This yields 
a direct measurement of electron density, independent of the state of 
equilibrium of the plasma. Klein [7] applied this method to two-dimen- 
sional flow around a model in an electromagnetically-driven shock tube, 
but the flow entering the standing shock was highly non-uniform. Wong 
[8] has utilized this two wave length technique to study the relaxation 
zone behind the moving shock in argon. His work, using the same combus- 
tion driven tube as the present study, reinforces the conclusion drawn 
from Figs. 2 and 3 that flow entering the standing shock is quite uni- 


am op aie 


The present study is part of a long range program at Stanford 
dealing with the dynamics and properties of shock-heated plasmas. In 
particular, the present problem has direct application to ionization 
phenomena in the close environment of hypersonic vehicles, such as 


communications blackout. 
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TI. EQUILIBRIUM FLOW THEORY 


A. Introduction 


In order to estimate the stagnation region electron density, it 
is necessary to calculate the state properties of argon at the stagnation 
line. Beginning with the state of the gas in the expansion chamber prior 
to the arrival of the incident shock wave (Fig. 1, region 1), the one- 
dimensional state calculations are made across the incident and standing 
Shocks, assuming adiabatic flow and the attainment of local thermodynamic 
equilibrium after each shock. The equilibrium assumption is justified for 
region 2 by the work of Wong [8], as well as the interferograms in Figs. 
2 and 3 of the present work. The cemputation is completed by assuming 
that the gas flows isentropically and incompressibly from the standing 


shock to the stagnation line. 


B. Gas Properties Across Shock Waves 


Properties across both shock waves are calculated using the Rankine- 
Hugoniot equations for a partially ionized, monatomic gas in conjunction 
with the Saha equation. The validity of the Saha equation for the 
equilibrium argon plasma has been demonstrated by Wong [8] and Cole [12]. 
Results of a computer program obtained by Wong were used to obtain the 
state of the gas in region 2. The equations used across the standing 


shock follow. 


The degree of ionization, x , is defined by 


N 
ee a 
~N, +N z 
A Me 


where Ny ee | : and N. are number densities of atoms, ions, and 
A 
electrons in any region of flow. The Saha equation, applied to the 


atom-ion-electron gas, is 
za 2 
3 x/kr 

B,(T) phy e 


= gal ‘ 
op, (t)(2m u)>/2(«r)?!? 


(21 pu 
A 











where p and T are the pressure and temperature, Ay is Planck's 

constant, x is the ionization potential, and wu is the reduced mass 

of the ion-electron pair. In order to make use of the Saha equation, 

it is necessary to approximate BA and B , wie DarlLEcion funerious 
A 

of the atom and ion. Since the neutral atom of argon has about 60 


excited states near the ionization potential, the partition function is 


approximately 
_ 162,500 
B(T)=1+60e * 


where 162,500°K. is the ionization temperature. Since the temperatures 
involved are less than 15,000°K., it is evident that the second term can 
be neglected. The argon ion has a ground state of multiplicity four, and 
a low level excited state of multiplicity two. Neglecting higher levels, 
the partition function of the ion is 


2060 
T 


B (T) = 4 wees. 
At 


Bond [9] estimates that neglect of these higher levels and the second 
lonization of the atom will not introduce appreciable error below 
15 ,000°K. Using these approximations, the Saha equation becomes 
3 G2 100 “1/2 
1° ile 

~ 2060\ ~ 
(21 13/2 (ny 9/2 plore 2 


ph 
Cae 


The equation of state of an ionized monatomic gas in thermodynamic 


equilibrium is 


(1+x) pRT = oRT ; 


ie. 
i 


where R is the average gas constant per unit mass of the mixture, and 


the masses of the ion and atom are considered equal. 


The one-dimensional equations of conservation of mass, momentun, 


and energy for flow across a normal shock are unchanged by ionization: 





2 2 
+ = 
PoVo Po P33 3 P2 
iy 2 se 
= Se 
5 Ve > 5 ie + Je 


The velocities (v) are given in the laboratory system, i.e. a system 
fixed with respect to the standing shock. These equations may be 


combined to obtain 


-ifi _ Lk E - 
h, is & | (> Pe) Creag. 


If excitation energy is neglected, the specific enthalpy is given by 


where e, is the internal energy per unit mass due to translation, and 


e. is the specific ionization energy. They are given by 


R*¥ 
Wo 


(l+x)T == R(1l+x)T  , 


and 


es ee = xRO 
ik p 0 m_(N_/x) 4 


Where R* is the gas constant per mole, W is the atomic weight of 
argon, W is the average mixture molecular weight, @ is the ioniza- 
tion temperature, and my is the mass of the argon atom. Finally, the 
specific enthalpy is 


(l+x)RT + xR@ 


roy 


ee 


Substituting this in Eq. (II-1), the Rankine-Hugoniot equation becomes 








Y p 
= ~ : x.-x : (II-2) 
es ie 
Pr 1+x,, Ty, 


A third equation, obtained by applying the equation of state on each 


Side of the shock wave, is 





cE x p fe) 
- 7 — = = (11-3) 
2 3 2 5 
A fourth independent relation is found by eliminating ue from the 
continuity and momentum equations to obtain 
"| = Ge 
ae: Sua. 


The densities are eliminated from the preceding equation by substitution 


from Eq. (II-2) and the equation of state. The result is 


p | p 
(3-3 
ow; ea Sa | 2s a i) 
(3-1) + 2( PE) E 
Py al +X, Ty, 


To obtain this equation in terms of Mach number, the speed of sound in 





Eero 2, a, a») LS, required. .Ffor.a monaronie, ont Zzedscasc ., a, is 


expressed as 


a= [7 Ft, - 3 R(14x,, )T, 


Dividing Eq. (II-4) by as , and solving for the pressure ratio P/Pp , 
the result is 


; 1/2 
p (x,-x,) OM 
3 at 2 32 Soe me 
BB 3 + 5MS + 5] (1-M5)* * 15 (iex,) TF, , 


where M,, is the Mach number entering the standing shock. This 





equation, along with the Saha equation, and Eqs. (II-2) and (II-3), 
form a determinate system for the four unknowns P2 ; as ee hate 


3 


x, -. The one-dimensional conservation equations used in this derivation 
are valid when applied to equilibrium regions on each side of a thin 
normal shock. Since the thickness of the non-equilibrium region in 
Fig. 3 is about half of the shock stand-off distance, the assumption of 


one-dimensional flow through the shock may be questionable. 


C. Stagnation Line Gas Properties 


If the flow from region 3 to the stagnation line is assumed adia- 
batic and non-heat conducting, then the total enthalpy is constant along 
any equilibrium streamline. In particular, along the central stream- 


isine, the enthalpy is given by 


he = h, + 


Nol re 


ye 
3 


The integrated one-dimensional conservation equations used across the 
shock are inappropriate as the flow becomes two-dimensional approaching 
the stagnation line. But Euler's equation along the central stream 


tube can be written as 


vdvo+ op = 0 ) 
0 
where neglect of viscous stresses is justified by the relatively high 
densities, hence Reynolds numbers, involved. The assumptions of adia- 
batic, non-conducting, equilibrium, non-viscous flow imply the assump- 


tion of isentropic flow. Integrating the previous equation, we obtain 


Pv 
ovdv + Oar ees ooag 


If the flow is also assumed to be incompressible, then p=p,=const. ; 


and the result is simply Bernoulli's equation: 


P= Py + 5 Pas 
Say = ato es 





The assumption of incompressible flow is hardly tenable for the hot, 
lonized gas, even at the low Mach number involved. But, since Rutowski 


[4] estimates the density change to be small, Bernoulli's equation will 


be used in these calculations. 


With stagnation enthalpy and pressure known, the previously derived 
equations of specific enthalpy and state, and the Saha equation are 
applied at the stagnation point to solve for QP, » x » and qT. 


These equations are 


Dy 
i 
Jn 


+ + 
5 ak x_)RT. x RO ’ 


Ps 


a (Tx JRE , 





=n 
ene 182,100 [2 
p n> a: 
Sp 
ae 3/2 5/2 S0Goa. 

nC ee ee 3-3 2) 

Results of the iterative solution were used to calculate N. from the 
S 

relation 

S A 
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ITI. APPLICATION OF INTERFEROMETRY TO ELECTRON DENSITY MEASUREMENT 


A. The Mach-Zehnder Interferometer 


A schematic of the interferometer is shown in Fig. 9. Monochromatic 
light from the source is collimated before entering the half-silvered 


beam splitter 5S Here the light is split into two beams, which travel 


a0 
via the fully silvered mirrors M or M, to be recombined at the 
second beam splitter 55 - pines light 1s viocused iby the Jens L,, » and 


the camera is placed at the focal plane of the test section. 


If the mirrors were all perfectly parallel, and the optical path 
lengths of the two beans exactly equal, one would expect reinforcement 
of the in-phase beams, resulting in a uniformly bright image at the 
camera. Now, if one of the mirrors is rotated slightly, a uniformly 
varying difference in optical path is introduced, resulting in aliter- 
nating reinforcement and destruction, i.e., a fringe pattern. If now a 
plasma is formed in the test section, the index of refraction, and hence 
the optical path length, of one beam is changed, resulting in a fringe 
shift given by 

A a male Ng 


Se ; 


dX dX 


where A (n-1) is the change in refractivity, L is the test section 
width, and A is the monochromatic observing wave length. The com- 
pensating chamber is introduced to compensate for the windows of the 


test section of the shock tube. 


Since the primary quantity desired is electron density at the 
stagnation line, the problem is to relate the observed fringe shift 
(from ahead of the standing shock, 2, to the stagnation line, s ) to 
the change in electron density. Conditions behind the moving shock, 2 , 
are taken from calculations by Wong [8], which were verified experi- 


mentally by him. 
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joe Development of Fringe Shift Equations 


In order to relate the observed fringe shifts to electron density, 
the refractivities of the neutrals, ions, and electrons must be related 


to thermodynamic quantities. 


The refractivity of the neutrals is obtained from the polarizabi- 


lity of the atom, a@, , by the equation 





‘A 
eT @ 
I, ee wile Nee ml, (Nam) = KymNy 


where K, = (21 a, /m,) is the Gladstone-Dale constant, which varies 


only slightly with wavelength. The atom number density is related to the 


degree of ionization by 


ie 


A (lax). 


=|? 


Tnerefore, the change in refractivity of the neutrals between any two 


regions is 
A (ny=2) = KA anes: tay me, 


Similarly, the ion refractivity is given by 


n l=2ra N can ae. el zs m.N 
PE tas ee aay a. 7 OG 
A is meh yee 2S, Cid GA 


The number density of the ions is 


OX 
N = N _ e 
+ 
A ee 


ot 
A(n ,-1) = et a5) ; 
at By 


10 


The refractivity of the electrons is given by Wong [8] to be 





we 
Al lg 
areal) ys 3 ’ 
Ww 
where Wp » the plasma frequency, is given by 
2 
wr = ea _ it a (ox) 
Pei oa mm p ? 
e eA 


w is the observing light frequency, and where the assumption has been 
made that the electron-atom collision frequency is much less than the 
observing frequency. This assumption is valid for conditions of the 
present study. The change in refractivity due to electrons is 
oe 
2 8 
A (n,-1) = S45 A (ox) = - 6.77 x 10°2"A (ox), 
om m7 
eA 
where » is in centimeters, and A (ox) signifies the difference in 


the product of density and degree of ionization between any two regions. 


Assuming that the refractivity changes for the various particles 


are independent, the total change in refractivity is 


AS = os (ny-1) +4 ns - 1) + A (n,-2) 


Recalling that the total fringe shift across, for example, the standing 


shock is given by 


L 
S = 5 A (n-1) j 


this fringe shift is expressed in terms of thermodynamic quantities by 


os “at 8 
A 
S=L¢5- pg(t-x,)-en{1-x,)+ g— (0.3% 3-P Xp) Onn OC lO X(6 2X 2"P Xp 


(Galea | 


ia 





The ratio of polarizabilities of the argon ion and atom is taken as 


oe 
A 
a = 0.67 , 


A 
and the test section width for the present experiment was L=5.0 em. 
Inserting these values and combining terms, the total fringe shift for 
ie OO A is 


y 


S = 2-70 510 [ose -8.86(p.x, : 0x) ] F Gone? 


oe 


where R500 = £1592 em?/ am has been used. Similarly, for A'=4030 A , 


Ky1030 = .1609 em? / pm » and the fringe shift is 


, 4 
oP => Gos x 16 [P.3-P5-T-13(0 x, - PpX,) | : (III-3) 


Combining the last two equations, the electron density in region 3 of 


the argon plasma is 


HI 


f snes i335 x 107" (s) = 2275) ; (IJI-4) 


3 2 


N 


where A= 4500 A, A’ = 4030 A, and S and S' are the measured 


fringe shifts at these wave lengths. With N, now determined, either 


5 
of the fringe shift equations, (III-2) or (III-3), can be used to solve 


for p Choosing the equation for X = 4500 A , the result is 


3° 


7 -5 ~22 
02 =P, 5.63 X10 ~ 1S 5.06% 10 (s = *., 


The degree of ionization can now be determined from 





The last three equations are sufficient to determine electron density 


and degree of ionization at any point in the plasma from observed fringe 


ie 





shifts between a known region and the region of interest. The fringe 
shift equations (III-2) and (III-3) may be used to predict fringe shifts 


from calculated plasma conditions, as in Fig. 4. 


13 
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IV. DESCRIPTION OF EXPERIMENT 


The experiment was performed using the Stanford aerophysics 
laboratory shock tube. Figure 8 is a schematic of the experimental 
arrangement. The 25 foot long driven section is constructed of extruded 
aluminum, with a two inch square inside cross section. The combustion 
driver uses a stoichiometric mixture of hydrogen and oxygen diluted with 
helium, which is ignited by a battery of spark plugs. The driver and 
driven sections are separated by a scribed aluminum diaphragm. Pumping 


> 


of the driven section, to about 10° mm. Hg before filling, is accom- 
plished through a station at each end of the section. The test section, 
located 20 feet from the diaphragm, supports the two-dimensional cylinder- 
plate model. A pair of optical windows permits viewing of flow about the 
model, and three barium titanate piezo-electric transducers just upstream 
of the model are used to measure the incident shock speed. The dump 

tank at the downstream end is separated from the driven section by a thin 


aluminum diaphragm which breaks on arrival of the incident shock. 


The Mach-Zehnder interferometer (Fig. 9), located at the test 
section, is on loan through the courtesy of Lockheed Missiles and Space 
Company. The light source was obtained by discharging 16 kilovolts 
from a 7.5 microfarad capacitor across magnesium electrodes spaced about 
1/4 inch apart. Although no direct measurement of the spark duration 
was made, oscilloscope pick-up of the capacitor "ringing" indicated 
that the spark was probably near maximum intensity during the entire 
exposure time of five microseconds. A blast shutter and pinhole were 
originally placed at the focal point of the light source in an attempt 
to reduce film exposure by plasma radiation. However, Fig. 1O shows 
that, even at relatively low Mach number, radiation reaching the film 
was intense enough to wash out completely the fringes through the 


5890A filter. 


Since it had become evident that the mechanical blast shutter was 
not fast enough for this application, a Kerr cell was obtained, again 
through the courtesy of Lockheed. This electronic shutter, manufactured 


by the Kappa Scientific Corporation of Sierra Madre, California, was 


14 





equipped with a pulse-forming network that permitted a very reproducable 
shutter-open time of five microseconds. For the present steady flow 
application, this exposure time proved satisfactory. Unfortunately for 
this application, it was designed for use in the ultraviolet, with 
acetronitrile as the working fluid. The Kerr cell was triggered to open 
for five microseconds during peak spark intensity, and the shutter-open 
period was always between 12 and 17 microseconds after passage of the 


incident shock by the model. 


The flow field was split by interference filters passing 4500 and 
5890 angstroms located just ahead of the Polaroid camera. Simultaneous 
measurements of near stagnation-line fringe shifts through the two fil- 
ters were sufficient to determine the electron density. However, 
variation in intensity of the dark fringes passed by the 5890A filter 
(see Figs. 2,3) indicated the presence of a harmonic. Further, compa- 
rison of observed fringe shifts with Eqs. (III-2) and (III-3) showed 
that the "5890" fringes actually corresponded to a wave length lower 
than 4500A . Examination of a spectrum obtained by Wong through this 
second order 5890A filter, also using a magnesium spark source, showed 
a third order pass band at 4030A . It was concluded that peaking of 
the ultraviolet Kerr cell at about 4000A , as well as the low intensity 
of the magnesium spectrum at 5890A , resulted in a fringe pattern domi- 
nated by ' = 4030A . The existence of harmonics in the fringe pictures 
indicates that the presence of some 5890A might cause a small increase 
in effective wave length above the value of A' = 4030A used in the 


calculations. 


An attempt was made to obtain total stagnation line heat transfer 
data using a platinum thick-film calorimeter, as described by Rutowski 
[4]. Initial efforts indicated that this phase of the program would 
require more effort than originally anticipated. Because of practical 
limitations of time, it was decided to concentrate remaining efforts on 


the interferometric studies. 


ats 





V. DISCUSSION OF RESULTS 


Usable fringe shift data were obtained for both A = 4500 A and 
X' = 4030 A on four runs. These were made at incident shock Mach 
numbers of 11.45, 11.8, 12.0, and 12.8. The lowest Mach number was 
obtained with an initial pressure (p,) of 10 mm. Hg. , while the others 
corresponded to Pj =2 mm. All tests were made in argon, with the a (ee 
inch diameter cylinder-plate model installed in the test section. Graphs 
of fringe shift, density ratio, electron density, and degree of ioniza- 
tion (Figs. 4,5,6 and 7) compare predicted and measured quantities in 
region 3, just downstream of the standing shock. Figure 7 also shows 
measured stagnation line degree of ionization and the corresponding 


prediction of equilibrium theory. 


Attainment of equilibrium in region 2, the flow entering the 
standing shock, was found by Wong [8] to occur about six microseconds 
after the incident shock (lab time) for the lowest Mach number of the 
present study, while the prediction of Petschek and Byron [10] is 
somewhat less. Since the present interferograms were taken between 12 
and 17 microseconds after the incident shock, it is reasonable that the 
flow entering the standing shock was in equilibrium. The interferograms 
in Figs. 2 and 3 show relatively straight fringes in region 2, supporting 
the conclusion of equilibrium there. The compensating effects of inci- 
dent shock wave attenuation and radiative cooling on the flow in region 
2 were discussed by Rutowski and Bershader [5]. Although the magnitude 
of these effects has not been determined, their effect probably was 
minimized by taking all the interferograms at the same time relative to 
the incident shock. Total available test time was estimated to be about 


100 microseconds for the fastest run. 


The interferograms at the two lower Mach numbers (Fig. 2a) show a 
sharp upward (positive) fringe shift across the standing shock due to 
change in refractivity of the heavy particles [Eq. (III-1)]. These 
shifts were found to correspond not to the predictions of equilibrium 
theory [Eg. (III-2)], but to the estimate obtained by assuming degree of 


ionization frozen at the value in region 2, and equilibration of the 
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translational modes only. Predicted quantities based on these assump- 
tions were obtained from the perfect gas shock tables of Keenan and 

Kaye [11], for y=5/3 , and the results are compared with measured 
values of fringe shift, density ratio, and degree of ionization in 

Figs. 4, 5, and 7. Since measurements at the stagnation line showed 
only a small increase in degree of ionization from that in region 3 
(Fig. 7), the flow field from the standing shock to the stagnation line 
is concluded to be definitely out of chemical equilibrium, and in fact 
nearly frozen. Measured values of electron density in region 3, which 
differed only moderately from those measured at the stagnation line, are 
shown in Fig. 6 to be much less than the corresponding equilibrium 


predictions. 


The interferograms at the two higher Mach numbers (Fig. 2b) show a 
radical change in the flow field. The shock stand-off distance is re- 
duced by more than a factor of two. The initial upward fringe shift 
again roughly corresponds to the frozen prediction, with the larger 
shift at the lower wave length, as would be expected from Eq. (III-1), 
for the heavy particles. But a large shift in the opposite direction, 
i.e. due to electron production [Eq. (III-1)] is also observed. Since 


this shift is proportional to wave length, it is greater for » = 4500A . 


The flow appears to approach local equilibrium at the line of 
maximum downward fringe shift, located about half way from the shock to 
ac Stagnation line. Since the fringe Shifty due to electrons Irom 
region 2 to 3 is somewhat less than predicted for equilibrium (Fig.4), 
the measured values of electron density obtained from Eq. (III-4) only 
approach the equilibrium estimate. It is possible that application of 
the one-dimensional conservation equations across a shock extending half 
way to the stagnation line is invalid, which could explain the lower 
measured equilibrium densities. A better theory would make use of a 
one-dimensional analysis across the frozen portion of the shock, and 
calculate the approach to equilibrium by a two-dimensional theory utili- 
zing a rate equation. However, the change in the fringe pattern (Fig. 3) 
perpendicular to the flow direction is small near the centerline of 


flow, indicating that the flow across the shock is nearly one-dimensional 
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in the region of interest. The flow from region 3 to the stagnation 
line causes an upward shifting of fringes, which for these conditions 
can only be caused by ion-electron recombination [Eq. (III-2)]. This 
effect, which of course is not predicted by the adiabatic equilibrium 
analysis, is illustrated by a small decrease in degree of ionization 
from region 3 to the stagnation line at M =1e.0 » and a larger decrease 
at the higher Mach number (Fig. 7). Although the Kerr cell effectively 
reduces the interferogram exposure from plasma radiation, other pictures, 
using only the blast shutter, show the existence of strong radiation in 
this Mach number range. It appears to be concentrated in the region 
behind the bow shock, and is probably largely due to radiative recombi- 
nation. It is concluded that the flow at the higher Mach numbers does 
approach equilibrium momentarily, but immediately begins cooling due to 
radiative losses. Radiative recombination could also explain the lower 
than expected equilibrium electron densities measured in region 3. The 
sharp change in stand-off distance over a very small Mach number range 
may be related to the attainment of equilibrium at the higher speeds. 


Table II summarizes these observations on the flow interferograms. 


Although no rigorous theory has yet been developed for convective 
stagnation region heat transfer in a partially ionized monatomic gas, 
the Fay-Riddell correlation equation [1] has been applied to this situ- 
ation by Reilly [3] and Rutowski and Bershader [5]. Predicted heat 
transfer rates based on this equation depend on the value chosen for 
the Lewis number, but do not vary markedly from the ideal gas prediction 
based on the same equation. Both Reilly and Rutowski and Bershader 
found substantial agreement of measured convective heat transfer rates 
with estimates based on the equation of Fay and Riddell. The present 
study suggests that the flow field ahead of Reilly's model was probably 
nearly frozen, since his study used a much smaller model and lower 
initial pressures than the present work, and involved Mach numbers of 
12 or less. A Similar comparison of the present study with the experi- 
ments of Rutowski and Bershader, involving Mach numbers greater than 
11.7, initial pressures of 10 mm. Hg. , and a larger model, indicates 


that equilibrium was probably attained in the stagnation region during 
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their work. These considerations suggest that the state of equilibrium 
of the shock-heated argon plasma does not critically affect stagnation 
region convective heat transfer. This result is not entirely unexpected, 
Since convective heat transfer is basically a boundary layer phenomenon. 
In the first approximation, it should be relatively insensitive to shock 
layer rate processes, particularly for the relatively high densities 


involved in the present study (Table I). 


Rutowski and Bershader suggest a radiative heat transfer model 
based on an extension of Kramers classical theory for free-free 
(bremsstrahlung) radiation to include radiative recombination into the 
upper bound states. Assuming thermodynamic equilibrium they estimate 
radiative heat transfer rate to be proportional to the square of elec- 
tron density. Use of predicted electron densities in the present case 
would have resulted in grossly overestimated radiative heat transfer 
for the frozen tests, but a much smaller error at the higher Mach 


numbers, where equilibrium was attained. 


The following factors were considered as possible contributors to 


error of the experimental measurements: 


le Error in determination of Mach number by the pressure 


transducer-oscilloscope system is estimated to be about 14%. 


oO. Actual initial temperatures varied from T =297°K 46 


z 
T=304°K . Since all theoretical calculations were based on 


T,=300°K » a small scatter in the (uncorrected) measured 


values would be expected. 


3. The work of Wong [8] indicates that the presence of impurities 
could cause a decrease in the relaxation time. However, this 
effect was predicted primarily for very low initial electron 
densities, where atom-atom collisions are predominant. Since 
all regions of interest in the present study involve relatively 
high electron densities, and the tube pressure before filling 

2 


was about 10° mm. Hg. for all runs, the effect of impurities 


is not considered significant. 
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Fringe shifts could be read to an accuracy of about one-tenth 
of a fringe. An error in this measurement would introduce 
only a small discrepancy between predicted and measured fringe 
shifts at a single known wave length (Fig. 4, A=4500A). 
Therefore, it is considered that qualitative observations from 
Fig. 4, i.e. separation of frozen and equilibrium regimes, are 


valid. 


However, the measured values of electron density, degree of 
ionization, and density ratio, which were all obtained from 
measured fringe shifts at two different wave lengths, are 
subject to more serious error. The proximity of the two ob- 
serving wave lengths (A=4500A, and A'=4030A) results in 
sensitivity of Eq. (III-4) to small errors in reading the 
fringe shift, especially at lower Mach numbers where the 
shift due to electrons is small. It will be recalled that 
original plans had called for using .A'=5890A as the second 
wave length, which would have resulted in considerably better 
accuracy. There is an additional uncertainty in the actual 
value of the lower wave length (A'=4030A + ~100A , - OA) 
This is derived from the possibility that presence of traces 
of 5890A might result in an "effective" A' > 4030A . 
Elimination of all 5890A radiation by procuring a suitable 
additional filter was precluded by time considerations. Use 
of a higher “effective wave length vice X'=4030A in the 
calculations would have resulted in a fairly uniform increase 
in the measured electron density for all runs. The estimated 
uncertainty in measured electron density due to these consi- 
derations is included in Fig. 6. It should be noted that this 
interferometric method is Capable of considerably greater 
accuracy than the present work obtained. However, the agree- 
ment of density ratios (Fig. 5) obtained by the two-wave 
length technique with the appropriate theory (frozen or equi- 
librium) corresponds to the conclusions drawn from the more 


accurate single-wave length fringe shift predictions (Fig. 4). 
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VI. CONCLUSIONS 


Flow entering the standing shock was in equilibrium during each 


GESt. 


Stagnation region equilibrium was not attained during the two 
lower Mach number tests. The degree of ionization remained 


essentially frozen at its value entering the standing shock. 


During the higher speed tests, local thermodynamic equilibrium was 
achieved ahead of the stagnation point. The thickness of the re- 
laxation region was about half of the shock stand-off distance. 
The degree of ionization after relaxation was less than the equi- 
librium prediction, and actually decreased from the end of the 
relaxation zone to the stagnation line. The lower than expected 
deésree of ionization in both regions 3 and s is attributed to 


radiative recombination. 


Since radiative heat transfer has been related to the square of 
electron density, the present work indicates somewhat less radia- 
tive heat transfer than predicted by equilibrium theory. But the 
observed recombination provides some justification for considering 
radiative recombination a primary component of stagnation region 


thermal transport. 
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TABLE I 
CALCULATED STATE PROPERTIES OF EQUILIBRIUM ARGON 
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X'=4030A 


A=4500A 





X'=4030A 


A=4500A 





Fig. 2b. Flow Interferogram (M_.=11.8) 
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A'=4030A 


A=4500A 


X'=4030A 


A=4500A 





Fig. 3b. Flow Interferogram (M_=12.5) 
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VS 
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Fig. 4. Fringe Shifts. 
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Fig. 6. Electron Density. 
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A=4500A 





Fig. 10. Blast Shutter Interferogram. 
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